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From these data the following points are to be noted: 

1. The vital index has the lowest mean value in the quarter ending on 
March 31, the winter quarter. In that period the birth incidence is rela- 
tively low and the death incidence relatively high. 

2. Next in value to this, but standing 9.40 ="= 1.75 points above it, is 
the mean vital index for the autumn quarter ending December 31. The 
difference being 5.4 times its probable error, may be regarded as significant. 

3. The spring quarter, ending June 30, shows the next higher mean, 
being 13.43 =*= 1.75 points above the winter quarter. 

4. The highest value of the index falls in the summer quarter, when 
births are most frequent and deaths least so. The mean value, however, 
lies only 5.54 =*= 2.05 above that for the spring quarter, a difference which 
cannot be regarded as significant. 

5. In variability of the vital index, the first two and the last quarters 
of the year, all exhibit significantly the same status. The vital index is 
distinctively more variable in the summer quarter, the difference in stand- 
ard deviations when this quarter is compared with that ending June 30, 
amounting to 6.04 ± 1.45. This may be regarded as significant, being 
4.2 times its probable error. 

It thus appears that the extremely close compensatory relation between 
birth rate and death rate, which Pearl and Burger 3 have shown to hold 
in annual figures, does not obtain within the single year. Instead there 
is a well-marked statistically significant intra-annual, or seasonal fluctua- 
tion of the birth-death ratio. 

1 Papers from the Department of Biometry and Vital Statistics, School of Hygiene 
and Public Health, Johns Hopkins University, No. 51. 

2 Cf. particularly Pearl, R., "The Vitality of the Peoples of America," Amer. J. Hyg., 
1, 1921 (592-674). 

3 Pearl, R. and Burger, M. H., "The Vital Index of the Population of England and 
Wales, 1838-1920," Proc. Nat. Acad. Scl, 8, 1922, pp. 71-76. 



A SOLUTION OF THE LINEAR MATRIX EQUATION BY DOUBLE 

MULTIPLICATION 

By I?rank Lauren Hitchcock 

DEPARTMENT OF MATHEMATICS, MASSACHUSETTS INSTITUTE OP TECHNOLOGY 

Communicated by E. B. Wilson, February 28, 1922 

1. Methods of Sylvester' and of Maclagan Wedderburn. — -A matrix A of 
order N is defined as a binary assemblage of iV 2 elements a ik , where either 
subscript may have any value from 1 to N. We add two matrices by add- 
ing corresponding elements. We multiply by the rule 

(AB) ik = Xa is b sk ; (s = 1, 2 A/). (1) 

If A u A 2 , . . . , A h and B{, B 2 , . . . , B h and C are all known matrices, 
while x is a required matrix, the linear matrix equation may be written 
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8(x) = A x xB x + A 2 xB 2 + . . . + A h xB h = C. (2) 

In 1884 Sylvester attacked this equation with characteristic energy. 1 
He pointed out that, since each element of the matrix 8{x) is a linear 
function of the elements of x — -the equation (2) being thus equivalent to 
TV 2 linear equations of ordinary algebra — it follows 8 has the properties 
of a matrix of order A r2 and satisfies a Hamilton-Cayley equation 

$ n - mrf"- 1 + m 2 e n ~ 2 - .. . + H) W = (3) 

where n = N 2 and where I denotes the matrix of order N 2 whose elements 
lik are unity when subscripts are equal, otherwise zero. If we know the 
coefficients nti . . . m„ we can solve (2), as is evident on multiplying (2) 
and (3) through by 0" 1 . We can find all the coefficients if we know the 
rule for finding m n which is the determinant of the matrix 8; for the left 
side of (3) may be found by forming the determinant of gl — 8 and after- 
wards putting 8 in place of g. To Sylvester, solution of the equation (2) 
was equivalent to finding, in terms of the matrices A\ . . . A h and B\ . . . B h 
the determinant of order N 2 belonging to the operation 8. 

His general method consists in setting up the two matrices of order 
N and of extent (Vendue) h 

x^A-i + X2A2 + . . . + x h A h , x x B x + x 2 B 2 + . . . + x h B h 
where x\, x 2 , etc. are scalar variables. By forming the determinants of 
these two matrices we have two quantics in the variables. He shows that 
the required determinant is a contrariant called the nivellant of these two 
quantics, and is an integral but not necessarily a rational function of the 
coefficients of the quantics. These coefficients are regarded as known 
quantities. The formation of the contrariant is in general difficult. 
He works out the solution for several special cases. With regard to the 
general case he remarks 2 

"Pour presenter l'expression geherale de ce determinant pour une matrice 
d'un ordre et d'une 6tendue quelconques, c'est-a-dire pour resoudre liqua- 
tion lineaire en matrices dans toute sa gen6ralit£, il faudrait avoir une 
connaissance des propri£tes des formes qui va beaucoup au dela. des limites 
des facultes humaines, telles qu'elles se sont manifestoes jusqu'au temps 
actuel, et qui, dans mon jugement, ne peut appartenir qu'a. l'intelligence 
supreme." 

But the indeterminate product of Gibbs was at that time unborn, or 
at least very young. 

In 1904 Maclagan Wedderburn 3 gave a solution of (2) by infinite series. 
His method consists essentially in selecting one of the known matrices, 
as Aj, multiplying by Aj 1 and into BJ 1 giving (2) in the form 

x + f(x) = A'jKCBj 1 = C 
where ^ is similar in form to 8 but contains h-\ terms. Then 
* = (1 + ^A)- 1 C = (1 - ^ + $ 2 . . . )C. 
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The choice of the term AjxBj to be freed from matrix coefficients affects 
the convergence of the solution. 

2. Transformation from Matrices to Double Dyadics.- — By Gibbs' double 
multiplication we may set up explicit formulas for the coefficients mi . . . m„ 
of the Hamilton-Cayley equation (3) . It is known that any matrix A 
is equivalent to a dyadic Saa' where a and a' are vectors. Any matrix 
term of the form AxB is equivalent to a sum of terms of the form aa'.X.bb' 
where X is the dyadic equivalent to x. We now introduce the double dot 
product defined by the rule ab:xy = a.xb.y, and may write 

aa'.X.bb' = ab'a'b:X. ' (4) 

Any term AxB is therefore equivalent to a sum of terms of the form 
MN:X where M and N are dyadics, and the linear matrix equation (2) 
is equivalent to 

?:X = 2MN:X = C (5) 

where <p is a dyadic whose antecedents and consequents are dyadics, say 
a double dyadic. 

3. Dyadics Treated Like Vectors in Space of Higher Dimensions. — -Adopting 
unit vectors e t . . . e M such that e;.e fe is unity when subscripts are equal, 
otherwise zero, any dyadic M may be expanded in terms of the dyads 
e;e fe ; and if we agree on some definite order, no matter what, among these 
dyads we may set E r = e,e fe , lettering r run from i to N 2 as i and k each run 
from i to N. With n =N 2 , any dyadic M ( may be expanded 

M ( = Eim/i + E2W,2 + . . . + E„m ln (6) 

where m a etc. are scalars. The double dyadic <p may be expanded as 

* = EiMi + E,M 2 + . . . + E M M re . (7) 

Now eie k :e r e s is unity when both i = r and k — s, otherwise zero. This 
is the same assaying E^tEj is unity when subscripts are equal, otherwise 
zero. It is immediately evident that the double dyadic <p will behave with 
reference to double dot product as would an ordinary dyadic in n- dimen- 
sional space with reference to ordinary dot product, and equally evident 
that the whole of the usual matrix theory will apply to <p with no essential 
change. 

4. Formation of the Hamilton-Cayley Equation. — Since <p:T/L = d(x) the 
Hamilton-Cayley equation for <p will be the same as that for 8. It re- 
mains to show how to calculate the coefficients wi . . m n of this equation. 
Ivet Mi, M 2 , . . . , M^, and Ni, N 2 , . . . , N^ be any dyadics: 

Definition.-— -The scalar ((MiNi, M 2 N 2 , . . . , M i ,N i ,)) s is defined to be 
the sum of terms computed as follows: the leading term is the product 

of double dot products Mi: N X M 2 : N 2 M^,: N^,; the other terms 

are of like form and are obtained from the leading term by keeping the 
antecedents M x .... M^ fixed in position while the consequents Ni . . . ~N P 
are permuted in all ways. The sign of any term is plus or minus according as 
the number of simple interchanges needed to form that term is even or odd. 



Vol. 8,1922 MATHEMATICS: F. L. HITCHCOCK 81 

Thus if p = 2 we shall have 

((MiWi, M 2 N 2 )) 5 = Mi:NiM,:N t - M,:N 2 M 2 :N!. 

If <pi, <p2, • ■ ■ , <Pp are double dyadics defined as in (7) the scalar ((pi, 
<pz, . . . , (pp))s is defined to be the result of expanding each <p as in (7), 
multiplying out term by term, and adding the scalars of all the terms of 
the product. 

If each M be expanded as in (6) the scalar elements of every <p will 
form a square matrix m ts of order iV 2 or n, so that each <p depends on N* 
scalar elements. The scalar ((^i, <&,.., <Pt))s is th e sum °f cubic determi- 
nants of order p whose main diagonals lie on the main diagonals of the 
respective matrices m ls ; the p matrices m ts are to be regarded as p non- 
signant layers of a cubic matrix. 4 

We may now suppose these p double dyadics to be all equal. The 
scalar just written becomes {(<p, <p, . . . to p factors)) s and may be abbrevi- 
ated ((/)) s . 

Theorem. — The scalar ((<p p ))s = p!m P where m p is the coefficient of 
ip n ~ t in the Hamilton-Cayley equation for <p. 

Proof. Let <p be expanded as in (7). The terms of the indeterminate 
product {<p, (p, . . . to p factors) are of the form 

(E^M,-, EjMj, E fe M„ . . . E r M r ) (8) 

where any subscript will have any value from 1 to n. Since E s : M ( = m ts 
the scalar of the expression (8) has for its leading term mumj^n^k ■ • • wi-rr 
and, by the definition, the scalar is the determinant of order p of which 
this is the leading term, a minor of the determinant of <p. Thus the scalars 
of all expressions (8) vanish except when the subscripts i, j, k, . . . , r are 
all different. Each choice of subscripts i, j, k, . . . , r will occur pi times. 
The scalar ((<p p ))s is therefore pi times the sum of minors of order p taken 
along the main diagonal of the determinant of <p. By the usual matrix 
theory then (0?*))s = p!m t as was to be proved. 6 

5. Transformation Back to Matrices. — It has thus been shown how the 
invariants mi ... w„ of the Hamilton-Cayley equation for <p or 6 may be 
calculated by products which are closely analogous with Gibbs' products 
of ordinary dyadics. To complete the solution we need to see how these 
invariants may be formed directly from the given matrices A\ . . . A h 
and Bi . . . B h without the need of first forming tp by the rule (4). • 

Consider first m x , the coefficient of <p n ~ l . If <p = SMN we have m t 
= SM:N = Sab':a'b = Sa.a'b.b' = 2A S B S where A. s and B s are the 
scalars of the dyadics corresponding to A and B, or, in matrix language, 
they are the respective sums of the elements of the main diagonals. Thus 

Mi = A 1S B 1S + A 2S B 2S + ... + A hS B hS . (9) 

For w 2 we have, by the distributive character of all the steps involved 
in the process, 
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2w 2 = ((*, <p)) s = S((M,N„ M»N*»s 
= -Z((8 i b' < a' J b ft a jfe b' )fe a'j fe bj fe )) s 

= -S[(a i b',:a' i b<) (a*b'»: a'*b») - (&ib\:a'» b k ) fob'*: a'A>») ] 
= Slaj.a'.-b^b'.-aft.a'fcbft.b'j: - a ! -.a' fe a fe .a' J b fe .b',b i .b' A ,] 

= 2[A iS Bi S A ks B kS - (^^ S ) S (B^) 5 ] (10) 

where, as before, A iS is the sum of elements of the main diagonal of A ir 
etc. and (AiA k ) s is the sum of elements of the main diagonal of the product 
of the two matrices A { and A k . 

By similar transformation we may show that 
6m 3 = •E[A is B iS A 1 - s B jS A kS B kS - A iS B iS (AjA k ) s (B k Bj) s 
- AjsBjsiAtAiUBiPds - AtsB^AiAjMBjBih 
+ (AiAjAJsiBtBjBds + (A,A^,) S (B^5,) S ]. (11) 

In every case all subscripts run from 1 to h, and terms where some or all 
of the subscripts are equal do not in general vanish. For a matrix term 
AxB is not in general equivalent to a single dyadic term MN : X. 

The general rule for any m is now easy to see. The quantity p!m p 
may be written as a summation of square bracket expressions, each square 
bracket enclosing p! terms. The leading term is 

2A iS B iS Aj S Bj S A kS B kS . . . A rS B rS 
where the p subscripts i, j,k, . . . r are those of the given matrix coefficients 
in (2) and may be either alike or different. The other terms are formed by- 
first keeping the antecedents A fixed and making interchanges among the- 
consequents B. A multiple interchange such that q consequents B a B b B c . . 
B e ByB g (whether alike or not) change places among themselves leads to a 
product of two factors as (A c A b A e . . . A d AfA a ) s (B a B f B d . . . B e B h B c ) s where 
the order of consequents is the order they stand after the interchange, while the 
order of antecedents is the reverse of that of the consequents. Any antecedents 
or consequents which have not been moved yield the same scalar as in the leading 
term. 

These quantities m may be transformed in many ways. 

6. A Second Method of Solution. — Another method of solving the equa- 
tion (2) follows at once from the analogy between 6 and an ordinary dyadic. 
If <p = SNM let <p' = SNM corresponding to 6', that is to *ZA'xB'.. 
Let C\, Cz, . . . C n „\ be n-i matrices or dyadics such that each is orthogonal 
to the given matrix C in the sense that C:C, = fori = 1, 2, . . . , n— 1. 
Form n — \ new matrices 6' (C { ) . A matrix orthogonal to all these new matrices 
will in general be a solution of (2) aside from a scalar factor which may be- 
found by direct substitution. This scalar may also be found as a function 
of m n and the matrices C, C u C 2 , . . . C„_i. 

This method is analogous with Hamilton's original method of solving 
the linear vector equation. 6 

7. Conclusion. — The analogy between a dyadic and a vector by virtue 
of the double dot product may be extended to polyadics and .fiT-tuple dot 
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products, hence to dyadics whose antecedents and consequents are poly- 
adics, say double polyadics. The analogy may also be applied to certain 
■other products. Thus the scalar ((ab, cd)) 5 or a.bc.d — a.dc.b may be 
regarded as a product of two dyadics M*N and dyadics found which are 
orthogonal with respect to star product, that is such that Mj*Nfc is unity 
when subscripts are equal, otherwise zero. We may similarly have star 
products for double dyadics and double polyadics. We shall have Hamil- 
ton-Cayley equations with respect to star products, where the scalars m 
are formed by star product just as those of this paper were formed by double 
dot product. The present paper is intended as introduction to more de- 
tailed treatment of some of these phases of multiple algebra. They have 
a bearing on ordinary algebra by virtue of the fact that a polynomial of 
degree K in N variables corresponds to a symmetrical polyadic of order 
K in space of N dimensions. 7 

1 In a series of Paris papers, C. R. Acad. Sci., 98, and 99 (1884) and elsewhere. The 
complete bibliography may be found in Bull. Inter. Assoc. Promoting Study Quaternions 
and Allied Branches, particularly for March, 1908. 

2 Paris, C. R. Acad of Sci., 99, p. 435. 

* "Note on the Linear Matrix Equation," Edinburgh Math. Soc. Proc, 22, 1904 (49-53). 
4 The laws of signant and non-signant indices for p-vr&y determinants were given by 
Lepine Hall Rice, Amer. J. Math., 40, No. 3, July, 1918. 
'By forming the determinant of *>-XI. 
'Elements of Quaternions, Vol. I, 2nd. Ed., p. 498. 

1 Binary quadrics as vectors in 3-space, with what I here call star orthogonality, have 
been studied in great detail by Emil Waelsch. See the bibliographies referred to in 
Note 1. His method is given in Wien. Ber., 112, 1903, p. 645, 1091, and 1533, and called 
"Binare-analyse." 



